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SUMMARY 

The flow field within the fuelslde preburner of the Space Shuttle Main 
Engine Is calculated using a reacting flow code (REAC12D). Inlet and modeling 
parameters Involved In the numerical calculation are systematically varied to 
establish the sensitivity of the calculated exit temperature profile. It Is 
found that differences In the Inlet equivalence ratio have a large effect on 
the turbine Inlet temperature profile. A variety of preburner Inlet modeling 
changes such as Inlet turbulence level, modeling the gases as burned, unburned, 
premixed, or unmixed, are shown to have a smaller effect on the calculated tur- 
bine Inlet temperature profile. Also, the form of finite differencing used is 
shown to have an effect on the temperature profile. 


INTRODUCTION 

The durability of the Space Shuttle Main Engine (SSME) can be a signifi- 
cant factor In establishing total launch costs of the Space Transportation 
System (STS), commonly referred to as the Space Shuttle. Frequent replacements 
of the engine or engine system components can potentially make the STS less 
cost effective than an expendable launcher. To reduce costs, NASA Is sponsor- 
ing an effort to Increase the durability of future SSME designs. This study 
Is part of that effort. 

The SSME design allows for a wider power range than previous rocket 
engines. This design calls for highly fuel rich or oxygen rich combustion In 
preburners to power the turbines that drive high pressure liquid oxygen or 
liquid hydrogen pumps and combustion is completed In a separate main combustion 
chamber. Due to flow losses and requirements of increased payload, the power 
level of the SSME has had to be Increased. This has caused more frequent 
Inspections and replacement of engines or engine components. One component 
that has experienced shorter life Is the fuelslde turbine which powers the high 
pressure liquid hydrogen pump. Durability of this turbine is strongly affected 
by the temperature profile leaving the fuelslde preburner. This temperature 
profile, though not accurately known. Is needed to design more efficient tur- 
bine blading and to make turbine blade life predictions. Experimental measure- 
ments of this temperature profile are difficult to make due to the severe 
operating environment within the SSME. Making measurements Inside the SSME 
would be expensive, time consuming, and could compromise the Integrity of the 
engine. Therefore, with the recent advances In computational capabilities, an 
easier design approach Is to use a reacting flow computer code to predict the 
Turbine Inlet Temperature Profile (T1TP). This was the objective of the 
current study. A two-dimensional code, REACT2D, was modified to analyze a 
preburner geometry with a hydrogen- oxygen single step reaction. A series of 
calculations were made to assess the sensitivity of the turbine inlet tempera- 
ture profile to variations in the flow entering the SSME preburner. 



SYMBOL LIST 


Cp \ specific heat for the 1th component, kJ/(kg - K) 

Cpm mean specific heat for the mixture, kJ/(kg - K) 

C v turbulence model constant 

f mixture fraction 

Hr heat of reaction, J/kg 

h enthalpy, J/kg 

k turbulent kinetic energy, m 2 /s 2 

M^ molar mass of 1th species, kg/mole 

mf mass fraction of fuel 

mi mass fraction of 1th species 

p pressure, N/m 2 

R universal gas constant, 8314.3 J/(kg - mol K) 

Rf rate of combustion of fuel, kg/(s - m 3 ) 
r radial distance, m 

S u source term for axial velocity, kg/(m 2 - s 2 ) 

S v source term for radial velocity, kg/(m 2 - s 2 ) 

S<p source term for the variable, (kg/(m 3 - s ) ) *( un Its of <p) 

T temperature, K 

u axial velocity, m/s 

v radial velocity, m/s 

x axial distance, m 

e turbulence energy dissipation rate, m 2 /s 2 

p density, kg/m 3 

v e ff effective viscosity: sum of laminar and turbulent viscosity. Ns/m 

Vfc laminar viscosity. Ns/m 2 

Pt turbulent viscosity. Ns/m 2 

f e ff the effective exchange coefficient for the variable <p. Ns/m 2 
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MATHEMATICAL FORMULATION 


The time averaged Navler Stokes (or Reynolds) equations for two- 
dimensional axlsymmetrlc flow are shown below (for details see ref. 1). 

Continuity: 

fc < r ' >u > * k (r <“ v > ■ 0 

Axial momentum: 

r [h < r ' >uu > * If < r '> vu >] - - « * f[lf ( r “eff li)* if ( r "eff if)]* s u 

Radial momentum: 

r [« <r ' ,uv) * If (r ' wv) ] * - f? * 7 [ii Kff if) * if ( r “eff if)] 

~ v ef f ~2 + S v 
r 

u cc the sum of the laminar and turbulent viscosities, 
ef t 

Scalar transport: 

f [if <w) * <ri>»*)] • f [if ( rr e ff i?)* if ( rr ef f i?) ] * \ 

r ff Is the effective exchange coefficient for the variable <p. 


The turbulent viscosity Is: 



Where k, the turbulent kinetic energy and e, the dissipation rate of kinetic 
energy are solved for using scalar equations. 

All but one of the reported calculations used hybrid differencing 
(ref. 1). Hybrid differencing Is a combination of central and upwind differ- 
encing. When the cell Peclet number (the ratio of convection coefficient to 
the diffusion coefficient) Is less than two, central differencing Is used. At 
higher Peclet numbers upwind differencing Is used. 

One calculation was done using Bounded Skew Upwind Differencing (BSUOS). 
While hybrid differencing uses a five point computation molecule, the more 
accurate BSUDS scheme employs a nine point computational molecule to more 
accurately calculate convection terms when the flow direction Is skewed with 
respect to the grid system (ref. 1). The scheme also uses a bounding routine 
to eliminate nonphysical oscillations or "wiggles" In the flow variables. 

The prechamber walls are modeled as adiabatic and law of the wall Is used. 
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COMBUSTION MODEL 


The stoichiometric reaction Is: 

H 2 + 2 °2 


H 2 0 f 1 .21x10® 


The enthalpy was defined as: 


h = m,H D + C T 
f R p,m 


The mean specific heat was given by: 

E"V 


P.m 


'Pi 


The specific heat for each component was approximated by a fourth order 
polynomial. The coefficients for the polynomial were calculated with a program 
by Mcbrlde and Gordon (ref. 2) using data from the Janaf tables, ref. 3. Two 
sets of coefficients are used for each species, one for a temperature range of 
100 to 1000 K, and another for 1000 to 3000 K to obtain better accuracy. These 
coefficients are given In table I. 

The eddy breakup combustion model of Magnussen and Hjertager (ref. 4) Is 
used. This model relates the combustion rate to the rate of dissipation of 
turbulent kinetic energy and concentration of species. Combustion Is modeled 
as a negative source term for the fuel mass fraction, mf. The Magnussen- 
Hjertager model was developed for both premixed and diffusion combustion, 
unlike previous models which were developed for only one type of combustion. 

In a diffusion type flame the oxidant and fuel exist In separate eddies, which 
must mix before combustion can occur. The time averaged reaction rate Is 
limited by the fuel concentration for lean flames, or by the oxygen concentra- 
tion for rich flames. In premixed flames the fuel and oxidant occur In the 
same eddies and the reaction rate Is determined by the flame spread, which Is 
related to the concentration of product. The Magnussen-Hjertager combustion 
model takes the limiting concentration of fuel, oxidant, or product and cal- 
culates a negative source term for the mass fraction of fuel. 


This source term Is: 
Sm 


P = -p ^ ml n Amp, 


m 


02 


STOIC 


AB 


m 1 

— 

S10IC ♦ ij 


A and B are empirical reaction constants and e/k Is the reciprocal 
of eddy lifetime. The present calculations used A = 32 and B = 4 unless 
noted otherwise. The quantities In the denominators are from the relations 
between product and reactants In the stoichiometric reaction. 

The rate of combustion Is then: 


R f - - S mf 


The combustion source terms were not calculated for Initial iterations to 
allow the computational flow field to be established. If this was not done, 
the solution could have diverged. As can be seen from the above equation for 
the combustion source term, no reaction will occur If the product concentration 
Is zero. The present program handled this by Introducing some product at a 
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specified Iteration In a specified cell or column of cells. Also, calculated 
combustion source terms could be so high that combustion was completed beyond 
what was physically possible. To correct this, the code was modified to check 
for and correct for these unrealistically large source terms. 

The species concentrations were calculated by using two scalar variables, 
the fuel mass fraction and the mixture fraction (f and mf). The mixture 
fraction was defined as the total mass fraction of fuel, both burned and 
unburned. 


f = m f + i8 m H ^ 0 

The program calculates the fuel mass fraction and mixture fraction each 
Iteration. Using these two variables other mass fractions may be found. 

The mass fraction of product Is: 

■y ■ 9(f - ra f> 

As there were only three species, the concentration of oxygen was: 

m 0 2 = 1,0 " m H 2 0 “ m f 

For this report, all the oxygen was consumed In the flowfleld unless a low 
value of A and B were used In the reaction model. 

The temperature may be found by manipulating the equation for enthalpy. 



The density was calculated using the equation of state for an Ideal gas 
mixture. 


GE0ME1RY 

The SSME configuration Is shown In figure 1. The fuel preburner Is 
located In the upper left side of the figure. The Inlet to the preburner, or 
preburner faceplate. Is three-dimensional. A section of the faceplate Is shown 
In figured (part of figure Is from ref. 5). The major portion of the flow Is 
through 264 oxygen and hydrogen coaxial Injectors. There are also a number of 
very small hydrogen Inlets used for area cooling. The computational grid used 
for these calculations Is shown In figure 3. Seventy axial and 53 radial grid- 
points were used. The preburner Is modeled as a series of coaxial nozzles 
Injecting Into a chamber. The walls are approximated by a rectangular mesh. 

The cooling flow which enters through the outer prechamber liner In the actual 
SSME was Included In the flow nearest the outer wall. The density of radial 
grldpolnts was Increased near the outer radius to ensure an adequate number of 
radial grldpolnts In the exit flow. This grid choice required CPU times of 
about 5 min on the Cray IS. The TITP was computed at an axial distance of 
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0.304 m, but the calculations were performed up to 0.5 m to Insure that exit 
boundary conditions did not affect the T1TP. 

The Inlet conditions were calculated from data provided by Rocketdyne. 
These are listed In table II. The Inlet was divided Into five zones by Rocket- 
dyne to calculate approximate temperature profiles. These zones were the 
Ignitor zone (In the center) and lettered zones corresponding to circular rows 
of Inlet nozzles. There are eight rows of nozzles corresponding to the let- 
ters A through H. Zone A- B Is next to the Igniter, followed by zones 
C-D, E-F-G, and H. Zone H Is the outermost zone and It Included the liner 
cooling flow. Inlet conditions used for the first calculation are listed In 
table III. Using a zero or very small Inlet velocity for the Inlet zone caused 
the program to diverge, so a nonzero Inlet condition was used for the Igniter. 
This flow was less than 0.6 percent of the total flow. Changing this flow 
slightly did not affect the TI1P. However If a kinetic reaction model were 
used, this might not be true. 

Originally, plans were made to run a vectorized three-dimensional code, 
after two-dimensional calculations, but the three-dimensional code has not been 
successfully vectorized. 

All of the two-dimensional calculations modeled the prechamber Inlet as 
coaxial Inlet zones. No walls were modeled In the Injector plane, lo model 
wall surfaces between Injectors would have Involved much more work In modifying 
the code. The prechamber Inlet has wall surface surrounding each Injector, as 
shown In figure 2. These surfaces cause a large velocity gradient which pro- 
motes higher turbulence. Higher turbulence could effect the calculated TUP. 

The outer . wall and prechamber dome were modeled as being adiabatic. All 
the calculations presented here required modifying these Inlet conditions or 
using combinations of them. 


RESULTS 

Baseline Calculation 

Using the Inlet conditions listed In table III, a baseline calculation was 
made. This calculation presents the closest approximation to the preburner 
design studied In this report. It Is provided as a baseline against which 
parametric variations can be compared. 

The TUP and velocity vector diagram for the baseline calculation are 
Shown In figures 4(b) and (c). The temperature profile shows much higher tem- 
peratures around r = 0.113 m. The TUP Is within the range of temperatures 
listed In table IV, which lists final temperatures for some mixtures. Also 
shown In figure 4(b) Is a TUP from a Rocketdyne calculation. The velocity 
vector diagram does not show any recirculation. Combustion Is very rapid for 
this calculation. The combustion causes the velocity to Increase by a factor 
of three within a few grid points which helps Increase the level of turbulence. 
Velocity only Increases again due to a reduction In flow area. 
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Modified Inlet Conditions 


The first modification was to change the inlet conditions to determine the 
effect on the TITP. The ignitor inlet condition for the baseline calculation 
was a mixture of combustion gases at a temperature of 1000 K. A second cal- 
culation with a cold (i.e., unreacted) Ignitor flow was made and the TITP was 
found to be similar. Another calculation which modeled the inlet gases as hot 
combustion products at the appropriate velocity and temperature was made. This 
calculation also gave the same TITP as the baseline calculation. However this 
calculation did show a small recirculation zone at the center of the preburner, 
on the prechamber dome (fig. 5). This series of calculations Illustrate that 
the inlet conditions assumed for the igniter can significantly affect the cal- 
culated flow field. Despite this effect, the TITP was insensitive to these 
changes. Apparently the ignitor represents such a small fraction of the entire 
flow that it does not significantly affect the exit temperature. 


Equivalence Ratio Changes 

A series of calculations in which the equivalence ratios of the inlet 
zones noted in figure 4(a) were varied. Changing the inlet conditions for the 
Ignitor zone did not effect the TITP. The long dashed line corresponding to 
changing the ignitor zone inlet condition to that of zone A B plotted on top 
of the TITP for the baseline (fig. 6(b)). The moderately long dashed line in 
figures 6(a) and (b) shows the affect of lowering the equivalence ratio of zone 
A-B and the Igniter to that of zone CD. The temperature profile was only 
moderately affected in the radial areas up to r = 0.112 m. Note, that the 
temperature is Increased by lowering the equivalence ratio because the equiva- 
lence ratio is so high. Changing the inlet condition of zone H to that of 
zone E-F-G only affected the outer portion of the TUP (short dashed line). 
However, changing the inlet condition of zone E-F-G to that of zone H had 
an effect over the whole TUP (see figs. 7(a) and (b)). Since inlet zone E-F-G 
accounts for almost half of the total flow this result should be expected. 

By using different inlet conditions TIlPs that were totally skewed in one 
direction could be generated. Using a set of inlet conditions with decreasing 
equivalence ratios resulted in a TI1P with temperature Increasing from the 
inner wall to the outer wall (fig. 8(b)). Keeping the overall mass flow simi- 
lar, a series of inlet conditions were found with increasing equivalence ratio, 
which yielded a TI1P that was skewed in the opposite direction. 


Diffusion Combustion Calculation 

All of the previous calculations used inlet conditions in which the oxygen 
and hydrogen were premixed for each zone. In the preburner, the flow is not 
premixed. The oxygen and hydrogen flow from separate nozzles. To model sepa- 
rate oxygen and hydrogen nozzles the number of inlet conditions was Increased 
(see table V). The mass flow rates of oxygen and hydrogen for each zone was 
kept the same as the baseline calculation. Also, the inner region of zone A B 
was approximated by a slower speed hydrogen nozzle to better model the actual 
geometry (fig. 2). The ignitor zone used the baseline inlet condition (mixed 
product and fuel at 1000 K). Temperatures vary widely downstream of the inlet 
nozzles due to the extreme variation in local equivalence ratio. For example, 
a temperature profile at x = 0.153 is shown in figure 9(a). Despite this the 
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TITP was very similar to the baseline calculation with only a minor variation 
In the outer radial portion of the TITP (fig. 9(b)). The temperature profile 
further downstream for this case came even closer to matching the baseline TUP. 
A color coded temperature plot for this calculation Is shown In fig. 9(c). This 
figure shows Intense mixing occurring near the rear wall of the preburner. 


The Reaction Model 

The reaction model uses turbulence quantities to calculate a combustion 
rate. To see If Inlet turbulence has an effect on the TITP two calculations 
were done using baseline Inlet conditions (five premixed zones). Increasing 
the Inlet turbulence for 5.5 to 30 percent Increased temperatures only slightly 
near the Inlet nozzles. Temperatures further downstream, and for the TITP, 
were unaffected. Decreasing the Inlet turbulence to 0.55 percent also had no 
effect on the TITP. Modifying the Inlet turbulence did not change turbulence 
levels to a corresponding degree throughout the flow field. Turbulence levels 
would better have been Increased throughout the flowfleld by larger differences 
In Inlet velocities, as was done for the diffusion calculation. 

As mentioned previously, the reaction model uses two empirical reaction 
constants. All of the previous calculations shown used A = 32 and 8=4. 
These numbers were quoted for a turbulent premixed combustion calculation. 

Using A = 4 and B = 0.5 which was quoted for a diffusion combustion calcula- 
tion did not produce much of a reaction In certain cases. Increasing A and 
B by a factor of two did not produce realistic temperature profiles either. 
Figures 10(a) and (b) show the effect of changing the reaction constants for 
baseline Inlet conditions and for Inlet conditions with no product. Decreasing 
A and B by a factor of two (from the high values) did not change the TUP 
for the baseline Inlet conditions but It produced a change In the calculation 
using no product In the Inlet. Further decreasing B for the baseline inlet 
conditions resulted In a lowering of temperature over most of the TIIP. Fur- 
ther decreases In the constants resulted In highly skewed profiles. Using 
still lower constants resulted In no combustion. 


Bounded Skew Upwind Differencing 

Finally a calculation was made using the more accurate BSUDS scheme. 

BSUDS is more computationally expensive and it Is more difficult to obtain a 
converged solution. To see If the BSUDS had an effect on mixing, a calculation 
was made modeling the Inlet as combusted gases. The TITP did differ from a 
hybrid differencing calculation as shown In figure 11. Whether or not this 
slight variation In temperature profile is significant to turbine blade life, 
the calculation Indicates the Importance of numerical accuracy in making react- 
ing flow calculations. The fact that a change was seen even with the large 
number of grldpolnts used In this calculation, provides strong evidence that 
fewer gridpolnts should not be used in these types of calculations. Routine 
three-dimensional calculations using such a fine mesh (using perhaps 70 by 55 
by 40) are currently Impractical from a time and cost standpoint. Less accurate 
three-dimensional calculations could be performed to establish gross trends in 
the flow field, but the results should be viewed by a knowledgeable designer 
who Is aware of the numerical Inaccuracies. 
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SUMMARY OF RESULTS 


The temperature profiles generated by these calculations with the eddy 
breakup combustion model show that the turbine Inlet temperature profile Is 
affected by changes In the equivalence ratio (If the affected area Is large 
enough), but that other aspects of modeling, such as modeling separate oxygen 
and hydrogen nozzles or modeling the Inlet conditions as hot combustion prod- 
ucts or unburned gases has little effect on the final exit temperature profile. 
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TABLE I. 


[C p = C] ♦ C 2 -T + C 3 -T 2 ♦ c 4 *t3 + C 5 -T 4 .] 



Ct 

c 2 

C3 

c 4 

c 5 

EMi 

EM 

4.06123 

3.57840 

3.29856 

-8 . 2329E-4 
-1 .1850E-3 
>1 .2650E-4 

2.9668E-6 
4 . 5000E-6 
3.2715E-6 

-1 . 2967E-9 
-3 . 34 54 E— 9 
-5 . 7020E-9 

5.5269E-14 
6 . 4591 E-l 3 
2.8912E-12 

T > 1000 K 

Mi 

In 

3.42384 

3.67138 

3.27290 

1 . 2083E-3 
5.4117E-4 
1 . 1 923E-4 

7.1516E-7 
4 . 5074E-8 
3.5771E-7 

-4.5067E-10 
-7 .9000E-1 1 
-1 .31 34E-10 

6 . 6681 E-l 4 
1 .5266E-14 
1 .4264E-14 


TABLE II. - ROCKETDYNE 
SUPPLIED INLET 
CONDITIONS 


Zone 

Mass flow, kg/s 

°2 

h 2 

A-B 

4.70 

4.81 

C-D 

8.46 

7.92 

EFG 

19.75 

17.44 

H 

7.52 

8.27 


TABLE III. - INLET CONDITIONS USED FOR BASELINE CALCULATION 
[An Inlet pressure of 250 atm was used.] 


Zone 

Mass flow, kg/s 

Radial distance, 
m 

Velocity temperature 

Equivalence 

ratio 

°2 

h 2 

m/s 

K 

IGN 

0.209 

0.2404 


77.6 

1000.0 

9.2 

A-B 

4.70 

4.81 


23.9 

162.8 

8.18 

C-D 

8.46 

7.92 


24.7 

162.5 

7.49 

EFG 

19.75 

17.44 


23.3 

162.4 


H 

7.52 

8.27 

.1088 - .1216 

25.6 

162.9 



































TABLE I V.v^ FINAL TEMPERATURE ASSUMING 
COMPLETE COMBUSTION 


Initial 

temperature 

Equivalence 

ratio 

Final 

temperature 

1000.0 

Infinity 

1000.0 

162.8 

8.18 

1076.2 

162.5 

7.49 

1144.5 

162.4 

7.07 

1191 .9 

162.9 

8.79 

1023.4 


TABLE V. - INLET CONDITIONS USED FOR DIFFUSION CALCULATION 


Zone 

Radial distance, 
m 

Velocity, 

m/s 

Mixture 

Number of 
cells 



- 


77.6 

h 2 * h 2 o 

6 


.0128 

_ 

.0192 

12.08 

h 2 

2 


.0192 

- 

.0224 

48.3 

h 2 

1 


.0224 

- 

.0288 

2.388 

°2 

2 


.0288 

- 

.0352 

48.30 

h 2 

2 


.0352 

- 

.0416 

2.388 

°2 

2 


.0416 

- 

.0448 

48.30 

h 2 

1 


.0448 

_ 

.0480 

46.95 

h 2 

1 



- 

.0544 

2.388 

0 2 

2 


.0544 

- 


46.95 

h 2 

2 



- 

.0672 

2.388 

°2 

2 


.0672 

- 


46.95 

h 2 

1 

EFG 


_ 

.0768 

44.26 

h 2 

2 


.0768 

- 


2.388 

0 2 

1 



- 

.0864 

44.26 

H2 

3 


.0864 

- 

.0928 

2.388 

0 2 

3 


.0928 

- 

.0992 

44.26 

»2 

4 


.0992 

- 


2.388 

0 2 

4 


.1056 

- 


44.26 

h 2 

2 

H 


_ 


48.98 

h 2 

2 


.1120 

- 

.1184 

2.123 

02 

4 


.1184 

- 

.1216 

48.98 

h 2 

2 
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Figure 1. - SSME powerhead component arrangement. 
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Figure 2. - Schematic of SSME fuel preburner. 
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Figure 5. - Velocity vector diagram for hot mixing gas calculation. 
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Figure 9. - Diffusion type combustion calculation. 
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Figure 10. - Different reaction constants. 
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